Density Perturbations and Black Hole Formation in Hybrid Inflation 
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We investigate the recently proposed hybrid inflation models with two stages of inflation. We 
show that quantum fluctuations at the time corresponding to the phase transition between the two 
inflationary stages can trigger the formation of a large number of inflating topological defects. In 
order to study density perturbations in these models we develop a new method to calculate density 
perturbations in a system of two scalar fields. We show that density perturbations in hybrid inflation 
models of the new type can be very large on the scale corresponding to the phase transition. The 
resulting density inhomogeneities lead to a copious production of black holes. This could be an 
argument against hybrid inflation models with two stages of inflation. However, we find a class of 
models where this problem can be easily avoided. The number of black holes produced in these 
models can be made extremely small, but in general it could be sufficiently large to have important 
cosmological and astrophysical implications. In particular, for certain values of parameters these 
black holes may constitute the dark matter in the universe. It is also possible to have hybrid models 
with two stages of inflation where the black hole production is not suppressed, but where the typical 
masses of the black holes are very small. Such models lead to a completely different thermal history 
of the universe, where post-inflationary reheating occurs via black hole evaporation. 
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I. INTRODUCTION 

A period of "inflation" or accelerated expansion in the 
early universe is an attractive idea in modern cosmol- 
ogy. Acceleration of the scale factor could drive the uni- 
verse towards homogeneity, isotropy and spatial flatness. 
However it is the ability of quantum fluctuations in the 
fields driving inflation to produce a nearly scale- invariant 
spectrum of quantum fluctuations that provides the most 
powerful test of the inflationary paradigm and may allow 
us to constrain the physics involved. Cosmological ob- 
servations allow us to measure the amplitude and tilt of 
the primordial density and, possibly, gravitational wave 
spectra on scales that would have left the horizon during 
inflation. 

The first inflationary models such as the old and the 
new inflationary universe scenario presumed that infla- 
tion began in the false vacuum state after the high tem- 
perature phase transitions in the early universe 
Later it was proposed that all possible initial conditions 
should be considered without necessarily assuming initial 
thermal equilibrium, and see whether some of these con- 
ditions may lead to inflation. This scenario was called 
chaotic inflation ||. For many years the idea of chaotic 
initial conditions seemed too radical, since it implied a 
considerable deviation from the idea of the hot Big Bang. 
It was argued that for a successful realization of infla- 
tionary theory one should satisfy so-called "thermal con- 
straints" jjj. However, gradually it was understood that 
the assumption of thermal initial conditions is neither 
natural nor helpful for inflationary theory |5j. As a re- 



sult, most of the models investigated now belong to the 
class of chaotic inflation, which provides the most general 
framework for the development of inflationary cosmology. 

The simplest models of chaotic inflation include theo- 
ries with potentials V((f>) such as m 2 </> 2 /2 or A0 4 /4. In- 
flation occurs in these theories at <j> > Mp. However, one 
may also consider chaotic inflation near cj> = in models 
with potentials which could be used for implementation 
of the new inflation scenario, such as —m 2 (j) 2 /2 + \<p A /A 
For brevity, one may call inflation in such models 
"new inflation" , to distinguish it from inflation at large 
4>, but strictly speaking these models also belong to the 
general class of chaotic inflation models: the original new 
inflationary universe scenario based on the theory of high 
temperature phase transitions have never been success- 
fully implemented in realistic theories. 

The simplest models of chaotic inflation such as the 



model 



/2 have many advantages, including natural 



initial conditions near the Planck density and the exis- 
tence of the regime of eternal self-reproduction of the 
universe Q. Normalizing the mass scale by the fluctua- 
tions in the microwave background observed by COBE Q] 
gives m ~ 2 x 10 13 GeV and the energy density at the end 
of inflation is V{4>) ~ (10 16 GeV) 4 . At this energy grav- 
itational waves contribute about 10% of the microwave 
background fluctuations. The tilt of the density pertur- 
bation spectrum in this model is n — 1 ~ —0.03. 

However, inflation occurs in such models only for 
4> > Mp. It is quite possible to have inflation at <j) > Mp 
in models with polynomial potentials, but in string the- 
ory and supergravity one often encounters potentials 
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which are extremely steep at (f> > Mp. It is not an un- 
solvable problem, see e.g. ||, but it would be nice to have 
a simple model where inflation may occur at <p < Mp as 
well. It is possible to achieve this, for instance, in ver- 
sions of "new inflation" with V{4>) ~ ~m 2 4> 2 /2 + A0 4 /4. 
However, in the simplest models of such type one has an 
unacceptably large negative tilt of the spectrum, unless 
the amplitude of spontaneous symmetry breaking is much 
greater than Mp J9|. Thus we return to the problem of 
having successful inflation at <p < Mp. 

There has recently been a lot of interest in the hybrid 
inflation scenario [|10|— |l5[ . Initial conditions for inflation 
in this scenario are not determined by thermal effects, 
and thus hybrid inflation also belongs to the general class 
of chaotic inflation models. However, hybrid inflation 
may occur at values of the scalar fields much smaller 
than Mp. The tilt of the spectrum in hybrid inflation 
typically is very small and positive, giving rise to so- 
called "blue spectra" Jl6|,0- 

The contribution of gravi- 
tational waves to the microwave background anisotropics 
is usually negligible. The reheating temperature in this 
scenario is typically large enough to ensure the possibility 
of electroweak baryogenesis, but small enough to avoid 
the problem of primordial gravitinos. 

It is still a challenge to obtain a natural implementa- 
tion of this scenario in the context of supergravity and 
string theory, but in globally supersymmetric theories 
this scenario appears in a very natural way. A very in- 
teresting version of the hybrid inflation scenario recently 
proposed by Randall, Soljacic and Guth was even called 
"supernatural" pl| . 

A distinctive feature of hybrid inflation is that it de- 
scribes the evolution of two scalar fields, (f> and ip. In the 
beginning one of these fields (field 4>) moves very slowly, 
and the second field may not move at all (though this sec- 
ond condition is not necessary [|l3|). The energy density 
supporting inflation is dominated by the false vacuum 
energy of the field ip. At the moment when the slowly 
moving field p reaches some critical value <p c , it triggers 
a rapid motion of the field ip, inducing a transition to a 
"waterfall" regime. Then the energy density of the field 
ip rapidly decreases, and inflation ends. 

Care is needed in evaluating the spectrum of density 
perturbations produced by inflation in the presence of 
more than one field. Many of the usual simplifying as- 
sumptions break down. Perturbations may no longer be 
purely adiabatic and hence curvature perturbations de- 
pend not only on the field fluctuations at horizon crossing 
but also upon their subsequent evolution up to the end of 
inflation, or even beyond. In the first versions of hybrid 
inflation the mass of the field ip whose false vacuum en- 
ergy density drives inflation was much larger than that 
of the slowly-rolling field <p an d so the single-field ap- 
proximation was quite sufficient. Also, inflation ended 
abruptly when the field p reached its critical value (j) c , 
and the field ip began its motion. 

This regime is certainly not the most general. Re- 
cently attention has been drawn to the possibility that 



both fields <p and ip could have masses close to the SUSY 
breaking scale m ~ 1 TeV, while the symmetry breaking 
scale in the ip direction may be as large as Mp [Q. In 
this case the process of rolling of the field ip towards its 
minimum may take a lot of time even if its mass is a 
few times greater than the Hubble constant at the end of 
inflation. This may be attractive in the context of super- 
symmetric theories but raises new issues about the gen- 
eration of density perturbations in the two-field model. 
We will consider these issues in this paper and spell out 
the dangers of ending hybrid inflation by a slow phase 
transition. 

The main problem associated with this scenario can be 
explained in the following way. The effective potential 
of the field ip used in |18) is symmetric with respect to 
the change ip — > —ip. As a result, at the moment of the 
phase transition the field ip can roll with equal probability 
towards its positive and negative values. This leads to the 
usual domain wall problem. To avoid this problem in the 
original hybrid inflation model [|l0| it was suggested to 
change the topology of the vacuum manifold and couple 
the field ip to gauge fields. In such a case instead of 
domain walls one may obtain either strings or monopoles, 
or (as in the electroweak theory) no stable topological 
defects at all. Monopoles should be avoided, but strings 
do not lead to any cosmological problems in theories with 
a relatively small scale of symmetry breaking, as studied 
in Ref. 10 1. Alternatively, one may consider versions 
of the hybrid inflation scenario considered in Ref. Jl3| , 
where no topological defects are produced. 

In the model proposed in Ref. Q topological defects 
do appear. In the simplest realization of this model 
one gets domain walls, which should be avoided at all 
costs. If one modifies the model to produce strings in- 
stead, one also has a problem, since strings corresponding 
to the scale of spontaneous symmetry breaking ~ Mp 
by themselves produce density perturbations Sp/p ~ 1 
on all scales. One could expect that monopoles would 
not lead to any trouble since the distance between them 
grows exponentially during the second stage of inflation. 
However, because of inflation, which occurs in this model 
during the long stage of rolling of the field ip to its min- 
imum, all topological defects in this model appear to be 
inflating, as in |l9|. Independently of the nature of these 
defects (domain walls, strings, monopoles, either topolog- 
ically stable or not) their exponential expansion leads to 
density perturbations 5p/p = 0(1) on the exponentially 
large scale corresponding to the moment of the phase 
transition. This may result in a copious black hole forma- 
tion. However, inflating topological defects in our model 
are rather specific, because they appear in the theory 
with \m^\ > H . For this reason the possibility of black 
hole formation by such defects requires separate investi- 
gation. This problem is extremely interesting since here 
for the first time the issue of inflating topological defects 
appears in the context of observational cosmology. 

Independently of this issue, the appearance of inflat- 
ing topological defects clearly demonstrates that the ex- 
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istence of the second stage of inflation in the hybrid sce- 
nario may lead to large density perturbations. As pointed 
out in Ref. p8| , the phase transition at <fr = cf) c leads 
to the appearance of a characteristic spike in the spec- 
trum of density perturbations. The existence of such a 
spike was first found in a similar model by Kofman and 
Pogosyan |20|. In the "supernatural" hybrid inflation 
model it is difficult to calculate the amplitude of the peak 
of the spectrum; in Ref. |l8| it was done slightly away 
from the point of the phase transition, where the ampli- 
tude of the density perturbations has already diminished. 
In order to perform the calculation, the evolution of the 
fluctuating field in was divided into several parts, and 
different approximations were used at every new step. 
However, the results of calculations of the amplitude of 
density perturbations near the narrow peak can be very 
sensitive to the choice of the approximation, especially 
in a situation where one may expect density perturba- 
tions to be large. Therefore we developed a more direct 
method of calculations of density perturbations in this 
model. 

One of the most interesting aspects of the model of 
Ref. |l8| is the existence of a regime in which quantum 
diffusion of the coarse-grained background fields domi- 
nates over its classical evolution and determines promi- 
nent features within our present cosmological horizon. 
Previously such phenomena were confined to scales much 
beyond our present horizon and were usually ignored. 
In Ref. Jig] , the machinery of stochastic inflation, see 
Ref. [^lj , was used to estimate the behavior of the fields 
close to the phase transition, where large quantum fluc- 
tuations make the stochastic formalism necessary. In this 
paper we will use this formalism to find whether or not 
most of inflationary trajectories come through the region 
where large density perturbations are generated. We be- 
lieve that the method which we developed may be of 
interest in its own right and can be applied to a more 
general class of models with many scalar fields. 

Our final results agree with the conclusion based on the 
topological defect analysis: density perturbations created 
at the moment of the phase transition are very large. In 
particular, in the model of Ref. |l^] with the parameters 
given there corresponding to the second stage of infla- 
tion lasting for 20 - 30 Hubble times one has den- 
sity perturbations Sp/p = 0(1). In such a situation one 
can expect copious production of huge black holes, which 
should lead to disastrous cosmological consequences. 

However, this is not an unsolvable problem. For a suit- 
able choice of parameters the second stage of inflation can 
be completely eliminated, and in this respect the model 
can be made very similar to the original hybrid inflation 
model of Ref. |Hj] , where the problem of black holes does 
not appear at all. A very interesting possibility appears if 
the second stage of inflation does exist, but is very short, 
lasting only two or three Hubble times H -1 . Then the 
black holes formed from the large density perturbations 
may be small enough to evaporate quickly. With the pa- 
rameters of the model of Ref. |18| the evaporation time 



is still very large even for the smallest black holes. How- 
ever, if one studies hybrid inflation with a larger Hubble 
constant, the black holes produced during inflation can 
be made very small, so that they evaporate before nucle- 
osynthesis. Even if the probability of formation of such 
black holes is suppressed, the fraction of matter in such 
micro black holes at the moment of their evaporation may 
be quite substantial, since the fraction of energy in radi- 
ation rapidly decreases in an expanding universe. This 
may lead to crucial modifications of the thermal history 
of the universe and may rejuvenate the possibility that 
the baryon asymmetry of the universe was produced in 



lere the second 



the process of black hole evaporation |2£ 

Finally, one may consider the models w : 
stage of inflation lasts for about ten Hubble times. As 
we will see, in this case the probability of the black hole 
formation may be sufficiently small, so that the amount of 
black holes does not contradict the cosmological bounds 
on the black hole abundance. It raises a very interesting 
possibility (see also |Q) that the black holes produced in 
the hybrid inflation scenario may serve as the dark matter 
candidates. In other words, dark matter may indeed be 
black! 

As we already mentioned, there exist some versions of 
the hybrid inflation scenario where topological defects do 
not appear at all. In this paper we will suggest another 
version of such a scenario, which we will call "natural" 
hybrid inflation. This scenario is a hybrid of the simplest 
version of "natural inflation" [^5 26 1, and the model of 
Ref. It shares some attractive features of "natural 

inflation" such as the natural origin of small parameters 
appearing in the theory. On the other hand, unlike the 
original "natural inflation" , our scenario does not require 
the radius of the "Mexican hat" potential to be greater 
than the Planck scale, which causes problems when one 
attempts to implement "natural inflation" in string the- 
ory p5| . We will show that in the models of natural 
hybrid inflation one can easily avoid the problem of large 
density perturbations. 

The plan of the paper is as follows. In Sect. 2 we will 
briefly describe the simplest hybrid inflation model (Tc| l 
and its relation to the model of Ref. pl|. We will find 
classical solutions describing the evolution of the fields 
<j) and ip in this model. Most of our investigation will 
be fairly general, but since the original hybrid inflation 
model jl0| is already well investigated [ ft2| ], we will con- 
centrate on the model of Ref. |18[ where an additional 
stage of inflation occurs after the phase transition. In 
Sect. 3 we will evaluate the amplitude of quantum fluc- 
tuations of each scalar field. In Sect. 4 we analyze the 
issue of inflating topological defects and the associated 
density perturbations. In Sect. 5 we study density per- 
turbations both before and after the phase transition and 
then discuss the important issue of quantum diffusion at 
the phase transition. In Sect. 6 we analyze the prob- 
ability of primordial black hole formation due to large 
density perturbations. We will also discuss the possibil- 
ity of reheating of the universe by evaporation of small 
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black holes. In Sect. 7 we propose and briefly describe 
the "natural" hybrid inflation model. We will discuss our 
results and summarize our conclusions in Section 8. 



II. CLASSICAL FIELD DYNAMICS 

The simplest realization of chaotic hybrid inflation is 
provided by the potential [frOf 



V{4>,^) = M 2 




+ irnV + l 7 0V- (1) 



For comparison, we will write here the effective potential 
of one of the models considered in Ref. [[181: 



V{(j>, ip) = M cos 2 



-fj) + I m 2^2 + I A 2^2^2 (2) 

y/2fj 2 v 4 V V K 1 



In the region ip < /,where inflation occurs in the model 
(g), the potentials (0) and (||) practically coincide, with 
the redefinition of parameters, f 2 — > M 2 /2\/A and A 2 — > 
27. (Note that the fields denoted and ^ in Ref. |l8| ] 
correspond to our fields "0 and 0. In this respect we have 
followed the notation of Ref. p2|.) In what follows we 
will study the model ([!]), but we will be most interested 
in values of the parameters close to those of Ref. |L8j . 

The equations of motion for the homogeneous fields 
arc then 



4> + iH4> = -(to 2 + ~fip 2 )<p, 

ip + 3Hip = (2\/AA/ 2 - 70 2 - \ip 2 )iP , 

subject to the Friedmann constraint, 



H 2 = 



8tt 
3M| 



f 



T/(0,^) + ^ 2 



+ of + 



(3) 
(4) 



(5) 



Although we can always integrate the equations of mo- 
tion numerically for arbitrary initial conditions, the clas- 
sical motion of the homogeneous field is not necessarily a 
good representation of the coarse-grained field on super- 
horizon scales. If the classical motion is sufficiently slow 
it can become dominated by quantum diffusion caused 
by wave- modes crossing outside the Hubble radius. 

In particular we wish to consider the case when cp is 
much larger than ip at early times so that ip has a large 
positive mass ~ jtp 2 and rolls rapidly to ip = 0.[] For 
values of <j> above a critical value <f> c , ip — is a stable 
minimum. Thus ip remains zero while <p slowly-rolls (for 



•C H 2 ) down to the critical value, 



where the 



ip field becomes massless. For smaller values, ip = is an 
unstable local maximum and quantum diffusion initiates 
a second-order phase transition from the false vacuum to 
the true vacuum state. 

In the simplest version of hybrid inflation, where the 
couplings A and 7 are of order unity, this is essentially all 
the dynamical evolution that matters. The bare mass of 



the ip field, 



-m 



2a/AM 2 , must be much larger than 



H jl2| and the phase transition occurs rapidly and infla- 
tion ends. The perturbation constraints on large scales 
are then readily derived from the usual single field results 
where the role of ip at early times can be neglected. 

But what if the bare mass of the ip field is not very 
much larger than HI In particular, what if this field 
has the type of potential we might expect for a moduli 
field with a minimum at ip ~ Mp and a negative mass- 
squared — to^ of order (ITeV) 2 about ip — 0? The false 
vacuum energy density at ip = is then M 4 ~ m 2 Mp ~ 
(10 n GeV) 4 and the Hubble constant H ~ ~ ITeV, 
as discussed in Ref. Q. It lies outside the range of 
parameters originally considered for hybrid inflation since 
it corresponds to an exceedingly flat potential for ip with 
effective coupling constant A ~ m^/M 4 ~ I0~ 30 . The 
roll-down of the ip field need no longer be fast and we 
must consider the complicated evolution of both the fields 
(p and ip during this stage. 

When ip = the potential simply reduces to V{4>) = 
M A + m 2 4> 2 /2. For the range of parameters and field 
values that we are interested in the constant term always 
dominates and the Hubble expansion can be taken to be 
de Sitter expansion with H = H = ^8ir/3M 2 /M P . It 
is useful to then write the bare masses of the two fields 
(p and ip relative to the Hubble scale as 



Hi 



and /3 = 2VA 



M 2 
H 2 



(6) 



In the case of a single scalar field a evolving during in- 
flation, one usually resorts to the slow-roll, a 2 <C V(a), 
and quasi-massless, V" H 2 , approximations to make 
analytic progress. This allows one to reduce the equa- 
tions of motion for the scalar field to a first-order equa- 
tion. However in our case the mass of the ip field is less 
than H only for a short interval, <p c \J\ — 1/p < <p < 
(pc^/TTTfp. Even then, we wish to consider values of a 
not much below unity so the quasi-massless approxima- 
tion may not be very good for <p either. 

Fortunately the fact that the potential energy, and 
hence the Hubble rate, are so nearly constant (and this 
really is a very good approximation for the parameters 
of Ref jl8|]) allows us to integrate the second order- 
equations in two regimes. 



*If instead <p rolls rapidly to <p = it remains there and 
the problem reduces to the usual case where the ip field rolls 
directly to the global minimum at V = 0. 
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Region I: small ip 



Region II: small <j> 



The first approximation regime will be for ip 2 / H 2 <C 
/3/7.Q This leaves the mass of the cp field constant and 
the equation of motion becomes 



cp" + 30' + a<f> = . 



(7) 



where a prime denotes a derivative with respect to N — 
Ho(t — t e ), the number of e- folds to the end of inflation. 
This can be readily integrated to give 



(j,(N) = exp(-r + 7V) + cp- exp(-r_ TV) . 



(8) 



The asymptotic solution is cp — cp + exp(— rN) where 
r = r+ > which approaches the slow-roll solution 
cp = cp + exp(— aN/3) for a«l. 

For cp > cp c , the ip field remains trapped in the stable 
minimum at ip = and we have effectively single-field 
inflation and the above solution gives the exact evolution 
of 4>(N). 

Below cp cp c we can no longer take the ip field to 
remain fixed at ip = 0. The ip field equation of motion is 



<//' + 3iP' = ( - 7-^5 ) , ■ . 



(9) 



We see that the field becomes massless at cp 2 = {(3/^)Hq. 
Re-writing this as an equation for ip as a function of cp 
gives 



dip 



^w +{1 ~ 2q)(t> Z + (kV +q2 ~ ^ = ■ (10) 

where 



2r 



- 3 - ^ 1 - 1 I^Vr 



r V 4 



(11) 



The exact solution is a linear combination of Bessel func- 
tions, 



V#) = <P q [-ciY v (k4>) + c 2 Mk4)\ 



(12) 



For <p <C (pc the growing mode is given by the small angle 
approximation 



V#) = -c x cP«Y y (Kq 



Cl 



Ci A< 



(13) 
(14) 



where the numerical coefficient A = (2 / k) v T \v) j 'it . 



As <p decreases, the effective potential for the ip field 
becomes dominated by its bare (tachyonic) mass vT^o- 
Thus for cp -C (p c , and still assuming ip 2 /Hq <C /3/A, wc 
have 



ip" + 3ip' - Pip = , 
which has the general solution 

iP(N) = ip+ exp(s+N) + ip- exp(s-N) , 



(15) 



(16) 



Matching to the asymptotic solution in Region I, 
Eq. (|l4|), we see that only the growing mode (s = s + > 0) 
exists in Region II. Recalling that N is measured from 
the end of inflation, it is simply given by 



iP(N) = ip e exp(sN) . 



(17) 



for all trajectories. The total amount of inflation after 
a given point is determined solely by the ratio ip e /ip, as 
noted in Ref. ]l^| . Even though the effective mass of ip 
has become large and negative, inflation will only end 
when ip 2 ~ V, which impliesn 



i> 2 * ipi - 



M 2 



4tts(s + 1) ' 
From a given value ip, the evolution will take 



s \sip 



(18) 



(19) 



e-folds to the end of inflation. As we shall see, if the 
initial value of ip is only of order Hq ~ ITeV there must 
be a further (32/s) e-folds before inflation ends. 

The growing value of ip increases the effective mass 
of the (p field through the interaction term and it also 
becomes larger than H . However, the cp field is much 
closer to the minimum of its potential than the ip field 
(4> C /Mp ~ 10~ 12 ) and <p soon starts to execute damped 
oscillations about cp — 0, as noted in Ref. Jig ], 

This time we have a parametric equation for the in- 
flaton field in terms of the triggering field, for which the 
general solution is 



<P(iP)=iP? [ciJ M (pV) + caY^piP)] , 



(20) 



where 



2s 



3 v^l 1 /9 r<m 

p = -—, p=— — , n=-\j--a. (21) 



s V 4 



^We also require ip 2 /Hg <C (3/X but in practice this is a 
much weaker condition for the parameters of |lq]. 



"'Note that for f3 > 1 this also ensures that de Sitter remains 
a good approximation until very near the end of inflation. 
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Matching to the asymptotic solution in Eq. (p. 3D, we 
find that only the J M solution is selected. We can see this 
from the small pip expansion of Eq. (EOh with C2 = 0, 



<P(ip)=c 1 iP p J^piP)^C!BiP 



(22) 



where B = (p/2) t2 /T(p + 1) and p + p = -(y - q) _1 = 
— r/s. Comparing this with the limiting behavior of 
Eq. (O) we find a relation between the coefficients 



ci B = (ci A) 



i/r 



(23) 



Therefore, as long as the solutions pass through the over- 
lapping region, our solutions of Region I evolve smoothly 
into Region II. 



III. QUANTUM FIELD FLUCTUATIONS 

In this section we discuss the evolution of quantum 
fluctuations of the fields. In the slow-roll approximation, 
the amplitude of quantum fluctuations of a massless field 
at horizon crossing (fc = aH) is approximately H/2n. 
However, in our case the masses of the fields are not 
necessarily much smaller than the Hubble scale, and cor- 
rections to the slow-roll result could be large. 

Since the potential V{(p,ip) ~ Vq to a very good ap- 
proximation, we can neglect the gravitational backreac- 
tion of the fields, and the equations of motion for linear 
perturbations in ip and <p can be written as 



Sip + 3HSip + — - (3H 2 + 7 2 Sip 



Sep + 3H5<P 



2-yip(PScp - 3Xip 2 Sip , 
k 2 



= —-fip 2 5(p — 2^ipcpSip . 



(24) 



(25) 



Note that in Region I, when ip = 0, the terms in the 
right-hand-sides are zero and the evolution of Sip and S(p 
decouple. We can write these equations, in terms of the 
canonically quantized fields u = aSip and v = aS(p, as 



k - I Ufe = . 



v'l + k 2 



Vk = 0, 



(26) 
(27) 



where primes denote derivatives with respect to confor- 
mal time, r\ = —1/aH, and we have chosen r/ = — 1 when 

cp = <p c . 

Since the mass of the <p field is constant, we can write 
an exact expression for the quantum fluctuations, 



2Vk 



^(-M^flW^-fe,), (28) 



where r = r + is defined in Eq. (||). This has the correct 
flat-space limit as — kr\ — > oo, Vk 
(p — > 0, and —krj—t 0, we find 



a ikr\ 



/v2fc, while as 



Vk(v) 



C(r) 



,i(l-r)7r/2 



2k 



(-kv) r 



where \U7\ 



C(r) = 2 



_ r r(3/2-r) 

r(3/2) 



(29) 



(30) 



This results in a scale-invariant spectrum of the growing- 
mode perturbations at horizon crossing with amplitude 



gfa = C(r) 



H 
2^ 



(31) 



Note that the coefficient C(r) gives a constant correction 
(independent of scale) to the usual amplitude of curva- 
ture perturbations (obtained in the slow-roll limit where 
C(0) = 1). 

On the other hand, the effective mass of ip changes 
with time and we cannot write down an exact solu- 
tion. However we can understand the qualitative be- 
havior by considering the effective Schrodinger equation 
for itfc, see Eq. (p6|), with the time-dependent potential 
v (v) = -(2 + (i)/rf + P/r] 2 ^ r) . This has a maximum 
value Vmax at r) max given by 



"max 



{fi + 2)r 
l-r 

13 + 2 
(3{l-r) 



/3(l-r) 



(3 + 2 

l/2r 



l/r 



(32) 
(33) 



On small scales, k 2 ^> V ma x, 
almost constant amplitude until 
it starts to diverge. As k — > oo we recover the constant- 
mass asymptotic (77 — > 0) solution 



the field Uk oscillates with 
-77 < y/(3 + 2/k, when 



C(- 



a i(l + S )7T/2 (_ 



2k 



(-k v y 



(34) 



where s — s + is defined in Eq. jl^). However, the am- 
plitude of the quantum fluctuations decay exponentially 
for k 2 < V(rf). Thus modes with k 2 < V niax will be sup- 
pressed. The amplitude of the growing mode of the field 
perturbations at horizon crossing can then be written as 



Sip* =C k ^- ■ 

Z7T 



(35) 



where the coefficient Ck is scale-dependent, as shown in 
Fig. 1 for the case (3 = 8. 

We can understand qualitatively the behavior of this 
growing amplitude of quantum fluctuations. For modes 
k ~ 1 that leave the horizon near <fi ~ <p c , the ip field 
is effectively massless and the amplitude of quantum 
fluctuations has the usual value H/2tt, with coefficient 



G 



Ck — 1. When < </> c (1 — l//?) 1 / 2 , corresponding to 
scales k > 1 + l/2r/3, the magnitude of the (imaginary) 
mass of the ip field becomes larger than H. Then the 
amplitude of quantum fluctuations even before horizon 
crossing, when k/a ~ m^, is Sip ~ m^/27r, which is 
already greater than H/2tt. As decreases and the cor- 
responding scale increases, the effective mass of the ^ 
field approaches its bare value while the amplitude of 
quantum fluctuations grows. At very large k we recover 
the constant-mass scale-invariant value, C^ — C(—s), 
see Eq. (|30|). For [3 = 8, this asymptotic value is 
Coo = 8.91 3> 1, see Fig. 1. On the other hand, for 



k < V m 



the mass (squared) of the ip field is large and 



positive at horizon crossing and the amplitude of these 
quantum fluctuations is suppressed. 



IV. INFLATING TOPOLOGICAL DEFECTS 

During inflation we must also consider the effect of 
short wavelength fluctuations that cross outside the hori- 
zon and perturb the coarse-grained background field on 
super-horizon scales. One can get a pretty good idea of 
the amplitude of perturbations on a scale corresponding 
to the time of the phase transition by investigating the 
inflating topological defects produced at that time. For 
this purpose one should study the evolution of the fluc- 
tuations Sip at the moment of the phase transition. 

Before the phase transition the field ip is very heavy 
and its quantum fluctuations can be neglected. As we 
have seen, its fluctuations are generated when the mass 
of the ip field becomes smaller than H , i.e. close to the 
phase transition. The exact duration of this stage is 
strongly model-dependent, but with the parameters used 
in Ref. |Q the time before the phase transition when 
the fluctuation can grow is only about if -1 . As a result, 
one may (approximately) visualize the scalar field ip at 
the moment of the phase transition as a sinusoidal wave 
with wavelength H^ 1 and amplitude Sip ~ H/^/2n, 



8ip c (x) 



H 

V2n 



sin Hx . 



(36) 



(Note that the amplitude here is H/y/2n, whereas the 
standard expression H/2it holds for the averaged ampli- 
tude (Sep 2 ) 1 / 2 , i.e. for the dispersion of the field.) 

This representation is not exact since many different 
waves give a contribution to Sip{x), but their wavelengths 
are comparable and, for the models we consider, the am- 
plitude of Sip in Eq. (|36|) is indeed of the order of H/ y/2ir. 
During inflation the wavelength of these perturbations 
grows exponentially, A ~ H a(t)/a(t c ), where a(t c ) is 
the scale factor at the moment of the phase transition, 
but the amplitude of the field S(p also grows exponen- 
tially. This is the main reason why spontaneous sym- 
metry breaking occurs despite the fact that formally the 
field ip, averaged over the whole universe, always remains 
equal to zero. During each interval of time At, when the 



effective potential can be represented as Vq — m 2 (t)ip 2 [2 
with m 2 (t) <C H 2 , the amplitude of the field Sip in ( |3q ) 



grows as 



Slpc 



H m 2 
■ exp ■ 



sin 



a(t)Hx 



whereas for m 2 (t) ^> H 2 one has 



At) 



Slpc 



H 

V2vr 



exp (m At) sin 



a(t)Hx 
a(t + At) 



(37) 



(38) 



In both cases spontaneous symmetry breaking during in- 
flation preserves the simple sinusoidal shape of the per- 
turbation, until the field rolls down to its minimum and 
a more complicated nonlinear regime begins. But this 
occurs already at the very end of inflation. 

Note that inflation near the cores of the topological 
defects continues for a while even after the field psi 
reaches its minimum at ip = / ~ Mp. Indeed, at 
that time the gradient energy density of the field dis- 
tribution @ ia - H 2 Ml exp(— 2N C ), where the fac- 
tor exp(— 2N C ) appears because of the stretching of the 
wave during inflation. For large N c it is much smaller 
than the potential energy density of the field ip, given by 
V(0) = 3H 2 Mp/8n. As for the kinetic energy density 
of the field ip it is always small near ip = 0. Therefore 
inflation in the vicinity of of topological defects contin- 
ues for a long time, until eventually the gradient energy 
becomes greater than V(0). 



As a result, after the field (38) approaches the mini- 



mum of the effective potential, on the exponentially large 
scale A ~ H _1 e we will have an extremely inhomoge- 
neous matter distribution. Roughly speaking, in half of 
the volume of the universe on that scale, the field will 
be near the minimum of its effective potential, whereas 
in the other half the scalar field will be close to the top 
of the effective potential with ip = 0, and inflation will 
be still going on. This shows that we will have density 
perturbations Sp/p = 0(1) on exponentially large scales 
corresponding to the time of the phase transition. 

This effect is very general. It is related to the inflating 
topological defects discovered in Ref. |[9| . The field may 
roll in any direction from ip = 0, but stable regions with 
ip = are constantly being created, corresponding to in- 
flating domain walls, strings, or monopoles. However, in 
Ref. jl9| the curvature of the effective potential was much 
smaller than H 2 . In our case the second stage of inflation 
occurs even if to 2 > H 2 , because the initial value of the 

field Sip c {x) = H/V2 n is much smaller than the ampli- 
tude of spontaneous symmetry breaking, so it takes a lot 
of time for the field ip to roll down. This specific feature 
implies that there will be no eternal self-reproduction and 
no fractal structure of topological defects in our model. 

Note that a similar effect may occur even in the models 
where topological defects are unstable if the decay rate 
of unstable defects is sufficiently small |2q| . A possible 
example may be provided by the metastable electroweak 
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cosmic strings or other "embedded defects" |29]. Such 
defects are unstable and usually do not cause any cosmo- 
logical problems. However, they also have ip = in their 
cores. If they inflate pq] and decay only after the end of 
inflation, they lead to large density perturbations in the 
same way as the topologically stable inflating defects. 

One may expect that large density inhomogeneities 
should lead to a copious black hole formation. The prob- 
lem of black hole formation in a hot universe is not trivial 
because pressure gradients do not allow black holes to be 
formed unless density perturbations on the horizon scale 
are of the order of one or even greater |j(],[3l],[l7 18 3^]33|] . 
But this is precisely the case for inflating topological de- 
fects. The centers of these defects remain forever at large 
density near the top of the effective potential, whereas far 
away from the cores the energy density gradually drops 
down to zero in an expanding universe. This corresponds 
to the delay of end of inflation SN 3> 1, and to density 
perturbations 5p/p > 1. If these topological defects were 
isolated objects immersed in empty Minkowski space, one 
can show that their thickness is much smaller than their 



Schwarzschild radius |19|. Therefore one would come to 
the conclusion that there is no suppression for the prob- 
ability of black hole formation in this scenario. However, 
it turns out that this is not what happens. 

To investigate this issue note that the sinusoidal wave 
( |36| ) is a particular example of a perturbation dip which 
is symmetric with respect to ip «-» —ip. The specific 
shape of these perturbations indicates that one cannot 
study them by considering Gaussian distributions cen- 
tered away from ip = 0. However, one could apply the 
theory of density perturbations to those parts of the de- 
fects away from ip = 0. The main ingredient is the cal- 
culation of the time delay for rolling of the field ip down 
to the minimum of its effective potential. For the centers 
of topological defects this time delay is infinite, which 
corresponds to density perturbations Sp/p 3> I. Thus, 
one may expect important black hole formation near the 
centers of inflating topological defects. However, as we 
will see, the mere fact that topological defects correspond 
to large perturbations does not automatically imply that 
they form black holes, because the distribution of density 
near the cores of these defects is extremely non-Gaussian. 
Note that topological defects correspond to places where 
the distribution of the field ip near their centers grows lin- 
early, ip(x) ~ H 2 x, see Eq. (|36j). Therefore the density 
distribution near the topological defects looks like a nar- 
row peak. Meanwhile, the distribution of the scalar field 
produced by the usual density perturbations at ip =/= 
has relatively flat regions where the spatial derivative of 
the field ip vanishes. These places have greater size, for 
the same amplitude of density contrast, and thus can be 
more easily converted into black holes. 

To estimate the density contrast Sp/p in the vicinity of 
the topological defects we will use the fact that during the 
last stages of inflation the field (p in our model is already 
close to the minimum of the effective potential. Thus, let 
us consider the field (B8I) at the moment when it reaches 



the minimum of its effective potential at ip = f. This 
happens approximately at a time t e given by equation 
(H/y/2n) e m ^ t,s = f. At this moment the distribution of 
the field ip, as a function of the distance r — xe~ Ht from 
the center of the topological defect, can be written in a 
very crude approximation as 



Sip c 



H 

V2V 



exp (m^, te) sm{Hr) . 



(39) 



Let us rewrite its evolution in terms of /, counting time 
from t e . For small r, 



Sip c (r) ~ fHr exp (m^ t) 



(40) 



which remains a good approximation until Sip c (r) ~ /. 
The time delay until the field at a distance r reaches / is 



SN{r) = HSt(r) 



H 



In Hr 



(41) 



This means that the density of the topological defect will 
exceed the average density by 5p/p ~ 5N(r) ~ 1 only in 
the core of the topological defect, at a distance 



rj-l ( m ip\ 

H cxp(-— J. 



(42) 



For the usual inflating topological defects considered in 
Ref. Jl9| one has m^/H -C 1, and therefore the devi- 
ation of density from the average density becomes large 
for the whole region from r ^ H~ x to r = 0. That is why 
one may expect such topological defects to look like huge 
black holes from the outside. Meanwhile in the model 
of Ref. Q one has m^/H > 1. As a result, the part 
of the volume inside the horizon where the energy den- 
sity deviates considerably from that on the horizon will 
be extremely small, being suppressed by exp(— Zm^/H). 
Therefore the total mass excess on the horizon scale is 
too small to cause the black hole formation. The same 
behavior of the density distribution is repeated at smaller 
scales as well. This suggests that inflating topological de- 
fects in the theories with m^p > H do not necessarily lead 
to dangerous black hole formation. 

The validity of this argument depends on many as- 
sumptions. For example, if the universe is matter dom- 
inated then black holes are formed even if the density 
perturbations are not very large. Therefore topological 
defects may serve as seeds for small black holes which 
are formed at soon after inflation, before the universe 
reheats and becomes radiation dominated. This may 
lead to very interesting consequences, to be discussed in 
Sect. 6. Meanwhile, black holes which are formed in the 
models with a prolonged second stage of inflation are very 
large, and they are formed very late, when the universe is 
supposed to be radiation dominated. Our investigation 
indicates that large black holes do not come from topo- 
logical defects. However, they may be produced by the 
usual density perturbations, which will be studied in the 
next section. 
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We should emphasize that our conclusions concerning 
black hole formation from inflating topological defects 
should be considered only as very preliminary. Indeed, 
here we deal with perturbations of density which are ex- 
tremely large, Sp/p ^> 1, and our simple analysis using 
methods of perturbation theory may be inapplicable. It 
is sufficient to remember that the interior of topological 
defects in our model continues to inflate for a long time 
after inflation ceased to exist in the outer space. Thus we 
have here a complicated wormhole-type geometry which 
requires a very careful analysis, similar to the numerical 
investigation performed in Ref. [p4| . It is the first time 
when we need to know a detailed structure of the universe 
filled by a gas of inflationary topological defects on the 
scale smaller than the size of the observable part of the 
universe. We need to return to investigation of this fas- 
cinating question in the future. In the meantime we can 
only say that it would be premature to conclude that the 
model of Ref. Jl8| contradicts observational data solely 
on the basis of our investigation of density perturbations 
produced by inflating topological defects. 

For this reason in the next section we will concentrate 
on the theory of usual (non-topological) density pertur- 
bations in models with two scalar fields and apply it to 
our model. The results that we have obtained from the 
investigation of topological defects suggest that, in the 
model under consideration, usual density perturbations 
may also happen to be very large. In what follows we 
will show that this is indeed the case. 



V. DENSITY PERTURBATIONS 

Quantum fluctuations of the scalar fields are responsi- 
ble for curvature perturbations on a comoving hypersur- 
face at the end of inflation, which can be evaluated as the 
change in the time (or number of e-folds) it takes to end 
inflation. In the case of adiabatic perturbations, e.g. in 
single-field inflation driven by the field ip, the amplitude 
of the curvature perturbation (on equal energy density 
hypersurfaces) remains fixed on super-horizon scales, so 
it can be calculated as SN = [H5ip/ip]* at horizon cross- 
ing where Sip can be estimated as H/2-k |3^|3f|. The 
origin of this curvature perturbation is that the jump of 
the field in the direction opposite to its motion leads to a 
time delay in the end of inflation which can be estimated 
by St — Sipjip. This equation immediately suggests that 
in the hybrid inflation model where the second stage of 
inflation after the phase transition begins at ip = 0, den- 
sity perturbations on the length scale corresponding to 
the moment of the phase transition should be extremely 
large, since at that time ip — 0. This is the standard sit- 
uation in all models where inflation occurs near the local 
maximum of the effective potential. It does not lead to 
any troubles if the corresponding scale is many orders of 
magnitude greater than the present size of the observable 
part of the universe. But in the model of Ref. (fj] this 



scale was supposed to be rather small, l c < 10 Mpc. One 
certainly does not want to have the observable part of 
the universe densely populated by huge black holes. 

However, these considerations are too naive. First of 
all, in our case we have two fields moving, <p and ip, so 
even if the field ip does not move at all, the whole field 
configuration evolves in time. Therefore the delay of the 
end of inflation is no longer given by the simple one-field 
expression St = Sipjip. Also, in the presence of more 
than one field, the perturbations may not be adiabatic 
and, as a consequence, their amplitude need not be con- 
stant [|7|j3£|. The effect of the perturbation must be 
integrated along the perturbed trajectory to the end of 
inflation, or until the evolution becomes adiabatic [ ]39"| , |38| | . 

Finally, because we no longer have a unique trajectory 
in field space, we must consider which trajectories are 
likely to be realized in practice. Because of the quantum 
fluctuations, the field ip does not stay exactly at the point 
ip = at the moment of the phase transition and we must 
investigate the dispersion of the probability distribution 
to see what the likely amplitude of density perturbations 
will be. 



A. Single-field perturbations 

Firstly we consider the simplest regime, before the 
phase transition, where single-field results apply. At large 
values of (p 3> 4> c the ip field has a large positive mass and 
remains fixed at ip = 0. The amplitude of ip fluctuations 
crossing outside the horizon is negligible. Thus we need 
only consider adiabatic fluctuations <50* along the trajec- 
tory, given by Eq. (31), which perturb the time it takes 
to end inflation. 



SN = 



H5<P 



C{r)H 
2'Kr6 Sr 



(43) 



The power spectrum of curvature perturbations on co- 
moving hypersurfaces, 1Z = SN , is then given by 



Tl 2 {N) 



g 0) 2 7 p 2r(N-N c ) 



(44) 



4tt 2 (3r 2 

where N c is the number of e-folds from the phase transi- 
tion ((p = <p c ) to the end of inflation (which we have seen 
could be of order 20-30). Assuming these perturbations 
are responsible for the observed temperature anisotropics 
in the microwave background, they give a constraint on 
the parameters of the model. The small (of order ITeV) 
Hubble constant during inflation implies that the con- 
tribution of gravitational waves to the microwave back- 
ground anisotropies will be negligible Jff],[f^]. The low 
multipoles of the angular power spectrum measured by 
COBE § give a value 1Z 2 ~ 3 x 10 9 on the scale of our 
current horizon, corresponding to 



Acmb — 46 



2 -m 

3 



M 



I0 1 



Gev) + 3 ln (l0 7 Gev) 



(45) 
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For parameters a and (3 or order one, this requires 
7 < 10~ 8 . Note that for these parameters we will have 
C > 10 4 TeV > H Q . 

This is one of the few cases in inflationary cosmol- 
ogy where we have an almost exact expression for the 
amplitude of curvature perturbations 40 1. The only ap- 
proximation we have made is to assume that the energy 
density remains constant, 2 /2 + m 2 <^> 2 /2 <C M 4 , so that 
we can neglect the backreaction on the metric. Using 
Eq. (E3|), we see that this will be true as long as 



8C(r) 2 M 4 
ZrK 2 ~Mj 



< 1. 



It follows that the allowed range of r is 



8 M 4 3 
3^ ^ r ^ 2 



8 M l 



9nTZ 2 Ml ' 



(46) 



(47) 



which for M <C 10 16 GeV essentially leaves a as a free 
parameter in the range < a < 9/4. At the upper 
limit, r ~ 3/2, the correction coefficient C(r) becomes 
large, giving a significant amplification of the curvature 
perturbations 1Z compared to the usual slow-roll approx- 
imation. 

Because the comoving scale at horizon crossing is just 
proportional to the scale factor oc e . the scale depen- 
dence of the power spectrum in Eq. (|4j) readily gives the 
tilt of the spectrum as 



1 



dlnTZ 2 
d\nk 



2r. 



(48) 



which becomes 2a/3 in the slow-roll limit. Note that in 
principle for this model one could have any value of the 
tilt in the range 1 < n < 4. A precise measurement of 
n, which may be possible with the next generation of 
satellite experiments EJ, would give a tight constraint 
on a. Present limits give n = 1.2 ±0.3 f?j|. Since n could 
be greater than one, there are also limits coming from 
production of black holes at small scales, n < 1.4 [ p~7| . 
Together they give a < 0.6, which constrains the size of 
the correction coefficient in Eq. (fji|) to lie in the range 
0.9 < C < 1. 



B. Two-field perturbations 

In this subsection we describe the interesting regime in 
which the system goes through the phase transition and 
quantum fluctuations of both fields become important. 
Here the curvature perturbation on a comoving hyper- 
surface at the end of inflation cannot be given simply in 
terms of that at horizon crossing since it does not remain 
constant on super-horizon scales. In order to compute 
the amplitude of the curvature perturbation at the end 
of inflation, one usually integrated the coupled differen- 
tial equations for the two fields' fluctuations and evalu- 
ated their amplitude at the end of inflation from that at 



horizon crossing. Only recently, Sasaki and Stewart |39| ] 
developed a formal method for computing the metric per- 
turbations at a given hypersurface from the change in the 
number of e-folds to that hypersurface as a local func- 
tion in field space, in the slow-roll approximation. This 
method was shown in Ref. (38| to be equivalent to the 
usual method of integrating the quantum field fluctua- 
tions. Unless one finds solutions for all trajectories in 
field space, the problem remains analytically intractable. 
In Ref. |3q | we found particular cases, with separable 
potentials for the interacting fields, in which the field 
trajectories were integrable and we could write explicit 
expressions for the amplitude of curvature perturbations 
at the end of inflation. In principle, for a general model 
all we need is a computer to evaluate the change in the 
number of e-folds to the end of inflation due to quan- 
tum fluctuations of the fields, for all points in field space. 
This gives us the possibility to investigate density per- 
turbations even for very complicated theories where it 
is not possible to express the final result in a compact 
analytical form. In the case of hybrid inflation we are 
fortunate to have a complete analytical solution for the 
classical evolution in Regions I and II which takes us 
from <f> = 4>c to the end of inflation. This will allow us 
to compute in a compact way the amplitude of curvature 
perturbations at the end of inflation (where the system 
becomes adiabatic and the amplitude remains constant 
on super-horizon scales) in terms of the field fluctuations 
at horizon crossing. 

The perturbation in the number of e-folds from any 
point (</>*,?/'*) in Region I to the end of inflation can 
be computed by evaluating the number of e-folds from 
a given point up until a surface of constant tp = ip m in 
Region II, since the time from tp m to the end of inflation 
is fixed, see Eq. (|l7j). Due to the overlap between the 
regions, this surface will also lie in Region I for ip m <C 
4>c = \J P/jHq. The number of e-folds to ip m is given by 



1 



where 

Cl, 



N(<t>*,<t> m )=N m - -In f , (49) 
l is a function of the trajectory parametrized by 



(ci Af'^x 



I s iU- r l s 



(50) 



Using the solution for tp(<p) m Region I, given by Eq. (p^), 
we can eliminate c\ to give 



s 



A 



(51) 



only in terms of the fields at horizon crossing. This is one 
of the few cases in which such an integration can be done 
completely up to the end of inflation, see also Ref. |38[ . 

Note that in general a trajectory beginning at a per- 
turbed point ((f) + S(p, ip + Sif)) may end up at a completely 
different point in field space compared with the nonper- 
turbed trajectory. This can make the comparison of the 
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lengths of the trajectories very complicated and could 
lead to entropy as well as curvature perturbations at the 
end of inflation. However, in our case all trajectories 
merge at the end of inflation, and this complication does 
not arise. The fact that by the end of inflation we are 
left with a single field (and thus all perturbations have 
become adiabatic) allows us to equate the amplitude of 
curvature perturbations on a comoving hypersurface at 
the end of inflation with perturbations on comoving hy- 
persurfaces at late times and, in particular, at the surface 
of last scattering, see Ref. Q . 

We can now evaluate the perturbation in the number 
of e-folds as 



SN 



l(k</>* 



S Y v (K(j>*) 



Sip* 
sip* 



(52) 



Note that Y v (z) ~ z~ v for small z, and thus the first 
term vanishes if the point (4>* , ip*) lies m Region II, which 
gives SN — Sip* /si})*, as required by Eq. (17). 

For <p* ~ cj) c we have dN/d(p* ~ l/(p c <C dN/dip*, for 
ip* <C 4> c - Since the quantum fluctuations 8(f>* and Sip* 
at this time are both of order H/2%, the amplitude of 
curvature perturbations is given by 



K = SN 



C k H 

2w sip* 



(53) 



where — 1, see Eq. (|35|). On the other hand, for 
(p* <C <fr c , the coefficient ft » 1, see Fig. 1, which is 
then responsible for large curvature perturbations. 

It is clear from Eq. (|53| ) that for small values of ip, the 
amplitude of curvature perturbations can become arbi- 
trarily large. In Fig. 2 we show a few equal-iV surfaces 
in field space (ip/ Hq, (p/(p c ), around and below <p = <p c , 
together with the line SN = 1. 

The amplitude of density perturbations on a comoving 
hypersurface when the curvature perturbations (^3|) re- 
enter the horizon is given by 



Sp 
P 



2 + 2w 
5 + 3w 



n. 



(54) 



where p = wp is the equation of state of the universe at 
re-entry. We therefore expect large density perturbations 
on scales associated with the phase transition. 



C. Quantum diffusion 

So far we have discussed the classical evolution of the 
homogeneous field and the effect of perturbations about 
the classical trajectories on a given scale for values of 
(p above (p c and below </> c . However at some stage the 
role of quantum diffusion of the coarse-grained field ip on 
super-horizon scales dominates over its classical motion. 
Purely classical trajectories in Region I beginning with 
ip ~ above <j> c are focussed along ip = 0, due to the 
large effective mass of ip at large </>, and continue to evolve 



close to ip = long after the point (p c when it becomes an 
unstable ridge. In practice we require quantum diffusion 
of the ip field to move the field off the ridge and begin its 
roll down to the global minimum. 

If diffusion washes out any trace of the classical motion 
as we cross <p = (p c , it does not make sense to calculate 
the curvature perturbations in terms of the classical tra- 
jectories. This would destroy our notion of associating 
points in field space with a given number of e-folds from 
the end of inflation. Quantum diffusion close to (p c could 
distort equal time hypersurfaces so much that we lose in- 
formation about the origin of trajectories. This problem 
is analogous to that of trying to trace the path of pho- 
tons beyond the (cosmological) last scattering surface. 
Beyond this surface, photons scatter many times and an 
observer at late times can no longer associate their path- 
lengths with a single smooth surface. Thermal diffusion 
is responsible for this loss of information in the trajecto- 
ries of photons beyond last scattering. In our case it is 
quantum diffusion of the scalar field that determines the 
loss of information. Note that in the above mentioned 
sense the region near the critical point <p = <p c becomes 
opaque to wavelengths equal to the wavelengths of per- 
turbations formed at <j> = <p c , but it remains transparent 
to perturbations with much greater and much smaller 
wavelengths. 

As argued above, our calculation of density perturba- 
tions relies only on being able to associate a given scale 
at late times, determined by the number of e-folds JV, 
with a unique smooth surface in field-space. Our re- 
sults for the amplitude of perturbations near <p = cp c 
in Eqs. ([32]) and (^) show that there is indeed a re- 
gion near <f> — (p c , ip — for which SN becomes very 
large. In Fig. 2 the region under the SN = 1 curve 
is the dangerous region. The question now is whether 
the parameters of the model are such that this affects 
a significant number of trajectories in field space. To 
determine this we should evaluate the probability distri- 
bution for the scalar fields and calculate how much of 
the initial wave-packet suffers large perturbations in the 
number of e-folds, i.e. SN > 1. A heuristic constraint 
could be that the size of the scattering region should be 
much smaller than the size of the wave-packet (very much 
like scattering of light of wavelength A by a target with 
a diameter a < A). These SN ~ 1 perturbations cor- 
respond to large curvature perturbations on asymptotic 
comoving hypersurfaces, which later become black holes. 
If a significant part of the packet enters the region where 
SN > 1, then at late times we cannot reconstruct the am- 
plitude of the initial perturbation, corresponding to large 
scales ((f) > (p c ). What happens is that, due to quantum 
diffusion, different scales will mix and their amplitudes 
will be undetermined for an asymptotic observer at late 
times. 

According to p^ ], in single field slow-roll inflation the 
regime of SN > 1 can be identified with quantum 
diffusion dominating over classical motion, 5<p > (p/H, 
i.e. with the well known self-reproduction regime p3j.p2[. 
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However, in two-field inflation this may no longer be the 
case. For example, fluctuations in one of the fields may 
not affect the time taken to end inflation. Even in our Re- 
gion II, where only ip determines N, it is the asymptotic 
time delay that determines SN, not the instantaneous 
perturbation [Sip/ip]» at horizon crossing. 

Let us now calculate the probability distribution for 
the field ip in Region I, as the field <f) = <f> c e~ r ( N ~ N °) 
slowly rolls down its potential. This can be done using 
the stochastic approach to inflation |2l],[l8| . Assuming a 
initial delta distribution for ip at <f> ^> <j) c , and an average 
quantum diffusion per Hubble volume per Hubble time 
rs i?/27r,^] the time-dependent probability distribution 
has the form 



1 



27TCT 



(55) 



where the dispersion a 2 it) satisfies the evolution equation 

Under a change of variables, x = exp[— 2r(N — N c )] and 
S(x) = a 2 (t)/H 2 , this equation becomes 



dS 



1 0(1 -x) 



dx 8ir 2 rx irx 
which has an exact solution 



S(x), 



1 / e 2 
o7C T \ ax 



(57) 



(58) 



where a = 0/3r ~ 0/a is a constant and T(a,z) is the 
Incomplete Gamma function. The solution S(N) char- 
acterizes the dispersion of the classical trajectories due 
to quantum fluctuations. Since the region where SN ~ 1 
has a width ip ~ H/2irs, see Eq. (^), at (j) = 4>c most 
classical trajectories will pass through this region. It is 
still possible that just one e-fold after the phase transi- 
tion the distribution will have spread so much that only 
a small fraction of all the trajectories still goes through 
this region. We thus consider the dispersion of the ip- 
distribution one e-fold after the phase transition, when 
x = cxp(— 2r). Note that the SN > 1 region is broader 
there, since Ck= e > 1 from Fig. 1, but the distribution 
( |58| ) is also wider. If the spread of the probability dis- 
tribution is still within the SN > 1 region at this stage, 



H 2 



exp(e 



-2r\ 



a e 



r(a,ae~ 2r ) < 



H 2 



An 2 s 2 



(59) 



§Note that our earlier analysis shows that this is an over- 
estimate for (j> > <j) c . However it should give a safe upper 
bound on the dispersion of the wave-packet. 



it will be difficult to avoid large perturbations on this 
scale. Note that this condition is totally independent of 
the coupling 7. We have plotted this condition in Fig. 3, 
where we show the contour plots of the dispersion of the 
distribution in units of the size of the SN = 1 region, 
(-0 2 ) 1 / 2 = nH/2-Ks, for n — 1, . . . , 7 in parameter space 
(a, 0). The figure shows is that in order for the distribu- 
tion to have spread (one e-fold after the phase transition) 
several times the size of the SN = 1 region one needs 
^> 1. The large separation of the lines indicates how 
difficult it is for the distribution to spread, unless ^> 1. 
If they were closely packed it would mean that, for not 
very large 0, the distribution would be much wider than 
this dangerous region and it would be possible for most 
trajectories to avoid this region. As it stands, for the val- 
ues of the parameters in Ref. 18 1, most of the trajectories 
will go through this region. 

In summary, we have shown that for a < 0.6 (the range 
of parameters allowed by observations of the spectral tilt 
on large and small scales) we find large perturbations, 
SN ~ 1, along most trajectories at the phase transition, 
unless /3> 1. This corresponds to large curvature per- 
turbations on these scales and thus to large density per- 
turbations after inflation which, as we shall see, leads 
inevitably to the formation of black holes. 



VI. BLACK HOLE PRODUCTION 

A. Probability of black hole formation 

We have seen that quantum fluctuations of the fields 
can be responsible for large curvature perturbations on 
a comoving hypersurface at the end of inflation. These 
perturbations, on scales that left the horizon 20 - 30 
e-folds before the end of inflation, as in the model we 
study, re-enter the horizon during the radiation era and 
could in principle collapse to form primordial black holes. 
The theory of production of primordial black holes from 
initial inhomogeneities was first discussed in Ref. f{0"|| , 
see also pl| . There is an expression for the probability 
that a region of mass m, with initial density contrast 
8(m) = Spj p\ m , becomes a primordial black hole, 



P(m) ~ 5(m) e 



-0 i /2S 2 



(60) 



where 2 ~ w p(i| |. In the derivation of this equation 
it was assumed that the universe was a barotropic fluid 
(p = wp) during gravitational collapse, and that the 
initial density contrast 5(m) is much smaller than one. 
There seems to be disagreement over the value of the pa- 
rameter and the way to calculate it in a radiation dom- 
inated universe, see Ref. |Q. Novikov et al. |32|,|33| give 
a simpler prescription based on numerical calculations. 
They claim that even for density perturbations less than 
one at horizon crossing, a black hole will form as long as 
the perturbation in the metric Sg a b = 2$ 5 a b is of order 
0.75 — 0.90, where $ is the gauge-invariant Newtonian 
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potential |33[. The analysis of a probability distribution 
for density perturbations with a peaked spectrum is be- 
yond the scope of this paper, but we expect a probability 
distribution like that of Eq. (|6p| ) with 1 % 1, somewhat 
larger than that of Rcf. 69]. Furthermore, Carr points 
out in Rcf. []30| that, for a scale invariant spectrum with 
density contrast 8 ~ 1 at horizon crossing, the probabil- 
ity of black hole formation is P ~ 1/2 on all scales, and 
half of the mass of the universe is always in the form of 
primordial black holes. From the above discussion it is 
clear that for the large density contrasts produced dur- 
ing the phase transition, 8 = (4/9) 8N ~ 0.5 at horizon 
crossing, see Eq. (|5^), there is no suppression of the prob- 
ability of black hole formation on scales associated with 
the phase transition. 

Note that in order to calculate the precise probability 
of black hole formation (|6^) , we have to solve the ambigu- 
ity in the value of the parameter (3. This would require a 
much more detailed investigation, see e.g. Q. However, 
for our purposes it is enough to realize that in a Gaussian 
distribution where the typical fluctuations of the density 
contrast are about 0.5, the fluctuations 8 ~ 1 are just 
one standard deviations away from 8 ~ 0.5. Therefore 
for every 10 horizon-sized regions with density contrast 
8 = 0.5 we typically find one region with 8 ~ 1, which 
can be expected to collapse and form a black hole. 

We should note that the previous discussion of the 
probability of the black hole formation is based on in- 
vestigation of conventional inflationary density pertur- 
bations. In our case, in addition to such perturbations, 
we have a dense gas of inflating topological defects. One 
might expect each of them to become a black hole, which 
would make the number of black holes much greater 
even than the one suggested by the estimates based on 
Eq. (|60|). (Here we are talking about the monopole- 
like inflating topological defects, since domain walls and 
strings with symmetry breaking / ~ Mp lead to a cos- 
mological disaster even if they do not inflate and form 
black holes.) However, we believe that this issue requires 
a more detailed analysis, see Sect. 4, and therefore in this 
paper we impose on the model of Ref. O] only those con- 
straints which follow from our investigation of the usual 
density perturbations. 

Let us now evaluate the typical size and mass of the 
black holes produced by these perturbations. Suppose 
that after the phase transition the universe inflated e N ° 
times. Then at the end of inflation the physical scale that 
left the horizon during the phase transition is Hq e Nc , 
where Hq, as before, is the Hubble constant during in- 
flation. Suppose that soon after inflation the equation 
of state became p = p/3, as for the ultrarelativistic gas. 
Then the scale factor of the universe after inflation grows 
as y/tHo. The scale H ~ 1 e Nc \ / tHo becomes comparable 
to the particle horizon ~ t at 



that time perturbations with density contrast 8 ~ 1 form 



t h 



(61) 



when the energy density becomes smaller than the infla- 
tionary energy density ~ H^Mp by a factor e^ 4Nc . At 



black holes of size H^e 2 ^ and mass 
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(62) 



For N c ~ 30 one would have black holes with mass ~ 
10 37 g, comparable to the masses of the black holes in 
the centers of galaxies. This is a very interesting mass 
scale, but copious production of such black holes would 
lead to catastrophic cosmological consequences. 

By changing the parameters of our model one can make 
the duration of the second stage of inflation rather short. 
For small black holes with N c = 0(1) in Eq. (|l]) should 
be modified because they are formed at the stage when 
the energy may still be dominated by the oscillations of 
the inflaton field with the equation of state p = 0. This 
changes a little our estimate for the time of formation of 
the black hole, 



and for the black hole mass, 



Af BH 
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(63) 



(64) 



The smallest black holes would correspond to N c ~ 1 
(and Ho ~ 10 3 GeV, as in |18|) would have a mass of 
about 10 11 g. Perturbations 8^1 giving rise to black 
holes in the mass interval 10 11 — 10 37 g are clearly ruled 
out by the bounds of Ref. jl7|. Thus, we should avoid at 
all costs the dangerous region SN ~ 1, since otherwise 
we will have too many large black holes. 

A very interesting possibility arises when one consid- 
ers such a peak in the spectrum, for not very massive 
black holes. From the bounds of Ref. |l7j we see that 
density contrasts of order 8 ~ 3/3/20 are just enough to 
give f?o = 1 in the mass range 10 15 — 10 30 g. Taking 
(3 ~ 1 from Ref. (3|] an d using Eq. (|65|), we find that 
a parameter s = 3 could indeed give the desired density 
contrast. This corresponds to a bare mass for the trig- 
gering field, m^p ~ 4 TeV, which is very natural. Further- 
more, the associated mass scale can be computed from 
N c = 32/3 - 11 as M B n = 2 x 10 20 g = 10~ 13 M Q . Us- 
ing the average density of our galaxy, p g ~ 10~ 25 g/cm 3 , 
we find that these small black holes may populate the 
halo of our galaxy and be separated from each other an 
average of 10 15 cm or about 6 times the size of the so- 
lar system. They could very well constitute the missing 
mass in our galaxy, and still pass undetected by the mi- 
crolensing surveys [Q . 

Note that changing slightly the parameters of the 
model one changes simultaneously the scale and the 
height of the peak in the black hole spectrum. This 
means that numerical values of the black hole masses 
and the distances between them can be made substan- 
tially different by modification of the hybrid inflation 
model. These numbers are very sensitive to the details 
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of the theory of black hole formation, which still requires 
a more complete analysis. It is important, however, that 
in the context of the hybrid inflation scenario the possi- 
bility that black holes may contribute to the dark matter 
of the universe becomes quite realistic. 

The idea that dark matter may consist of black holes 
produced after inflation was explored earlier by Ivanov, 
Naselsky and Novikov |24| . They performed a detailed 
investigation of the probability of formation of large black 
holes in such models, and in this respect their work can 
be extremely useful. Their model required the existence 
of a plateau in the effective potential, which would lead 
to a high peak in the spectrum of density perturbations. 
However, it is very difficult to obtain a realistic model 
of a single scalar field where one has almost exactly flat 
spectrum Sp/p ~ 5 x 10~ 5 on all scales from 10 28 to 10 20 
cm, and a sharp peak with Sp/p ~ 1 on a slightly smaller 
scale. Meanwhile, as we have shown, in the hybrid in- 
flation scenario this possibility emerges in a very natural 
way. 



B. Reheating from black hole evaporation 

A very interesting application of the above results 
comes when we consider a two-stage inflation with a suffi- 
ciently short period of inflation after the phase transition. 
With the parameters of the model p8| even the smallest 
black holes are very heavy and evaporate too late. How- 
ever, by choosing a model of hybrid inflation with a suf- 
ficiently large Hubble constant and short second stage of 
inflation one can have a very interesting regime when the 
black holes will dominate the energy density of the uni- 
verse soon after the end of inflation and later evaporate 
before nucleosynthesis, reheating the universe. 

Let us first assume that small black holes were formed 
in a radiation dominated universe, soon after the usual 
stage of reheating after inflation. To evaluate the proba- 
bility of black hole formation we need to know the density 
contrast at horizon crossing during the radiation domi- 
nated stage, see Eqs. (p3|) and (|54|), 
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The number of e-folds in the second stage of inflation 
is given by N c — (l/s)\n(ip e /ip). With initial condition 
ip ~ Hq/2it we have 



1/s 



(66) 



The time it takes a black hole to evaporate is given by 



1 



g*M P \ Mp 
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(67) 



Here g* ~ 10 2 is the effective number of particle species 
at the time of the black hole evaporation. For black holes 



with small N c formed at the matter dominated stage we 
have 



g*M P \H 



(68) 



Suppose that the fraction of matter in the black holes 
initially was only very small, and the universe was radi- 
ation dominated from the time of black hole formation 
to the time they evaporate. Then the fraction of mass in 
black holes grows during this time as a(t) ~ \fi due of the 
more rapid decrease of the energy density of relativistic 
particles outside black holes. At the instant before the 
black holes finally evaporate, the fraction of energy in 
black holes has grown by a factor of 




AT c /2 



(69) 



Therefore even if only a small fraction of energy was in 
the black holes initially, because the probability of their 
formation was suppressed by the exponential factor in 
Eq. (|6C|), we only require P(8) > (r jth)~ x l 2 for the black 
holes to come to dominate the energy density of the uni- 
verse before they evaporate. 

To give a particular example, let us consider hybrid 
inflation with Ho ~ 10 14 GeV (which is much greater 
than in the model of Ref. ps[ ). Let us take, e.g., s ~ 3, 
i.e. larger than the usual parameters of the models of 
Ref. [|l8[ but much smaller than those of Ref. (lO) . Then 
we have the total number of e-folds at the second stage 
of inflation e ° w 40. The density contrast in this 
case is S ~ 1/6. The fraction of matter in the black 
holes, according to (|60|), will be about 10~ 5 . In fact, 
this is a rather conservative estimate, since our investi- 
gation of topological defects suggests that this number 
may be much greater. The black holes will be produced 
at the moment th ~ 6 x 1CP 34 s. They will have mass 
Mbh ~ 6 x 10 4 g, and will evaporate at t ~ 2 x 10~ 16 s, 
much earlier than the epoch of nucleosynthesis. Because 
of the large growth of the scale factor and redshift of en- 
ergy of relativistic particles, at the time of the black hole 
evaporation practically all matter in the universe will be 
in black holes. This means that practically all particles 
which exist in the universe at t > 10~ 16 s are created at 
the moment of the black hole evaporation. 

Note that in the above example, although the prob- 
ability of black hole formation is very small, they still 
give the dominant contribution to the energy density at 
late times because the energy of relativistic particles de- 
creases much faster than that of black holes. However, 
this black hole dominance may begin much earlier if they 
are formed before conventional reheating is complete and 
the equation of state is p ~ 0. This condition can be eas- 
ily satisfied in the case of very small black holes. Then 
there will be no exponential suppression of the probabil- 
ity of the black hole production ( ]6C| ) , and the fraction of 
energy in the black holes could be large from the very 
beginning. 
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The process of black hole evaporation could be respon- 
sible for the baryon asymmetry in the universe, even 
though it is not very easy to get large baryon asymme- 
try by this mechanism |2q] . Typically it is assumed that 
reheating and thermalization of the universe occurs due 
to the inflaton field decay and the subsequent particle in- 
teractions, or through bubble collisions like in first order 
inflation. The natural assumption was that the gravita- 
tional interaction at the stage of reheating could be ne- 
glected. Here we have another, very unusual mechanism 
of reheating. Even in the absence of bubble wall collisions 
or a large coupling of the inflaton to matter, a consider- 
able fraction of matter after inflation could be in the form 
of small black holes. Unlike in the extended inflation sce- 
nario jiH), in our case all such black holes are formed in 
the same mass range given by Eqs. ([32|) and (|64|). If the 
probability of black hole formation is not strongly sup- 
pressed, then very soon they dominate the energy density 
of the universe |E3] . Eventually the evaporation of these 
black holes could reheat the universe. This opens up an 
interesting possibility of connecting the origin of matter 
in the universe with black hole physics. 

Let us estimate the reheating temperature of the uni- 
verse in this scenario. Black hole masses in the process of 

1/3 

their evaporation decrease as Mbh [l~V T ] • (Here we 
have taken into account that the age of the universe at 
the moment of the black hole formation is much smaller 
than their evaporation time r.) The main part of the 
energy release by the evaporating black holes occurs at 
the end of the time interval r. Therefore one may simply 
use the standard temperature-time relation for the hot 
universe to get the following estimate of the reheating 
temperature T r after the black hole evaporation: 

For a particular example which we studied (H ~ 10 14 
GeV, s ~ 3, g* ~ 10 2 ) we get T r ~ 2 x 10 9 GeV. This 
estimate is extremely sensitive to the choice of the pa- 
rameters. One can easily get reheating temperature as 
high as 10 10 GeV or even greater, or as small as 1 eV. 
The only real constraint on this temperature is that one 
should be able to produce the baryon asymmetry of the 
universe during or after black hole evaporation, and be- 
fore nucleosynthesis. This picture differs considerably 
from the standard theory of reheating due to the decay 
of the inflaton field, see e.g. |46|j47|] . 

Perhaps one can appreciate a potential importance of 
this regime if one remembers that the standard reheating 
due to the inflaton decay often is very inefficient because 
of the small coupling of the inflaton to matter |^,|| . In 
such cases the universe for a long time remains matter 
dominated (p = 0). In some other cases reheating is 
extremely efficient in the very beginning, but later be- 
comes inefficient, so that the universe eventually may 
become matter dominated again p7| . But then forma- 
tion of black holes is no longer suppressed by radiation 



pressure. In this case matter easily collapses into small 
black holes, which later evaporate and reheat the uni- 
verse. Thus, the absence of the usual reheating triggers 
black hole formation, which eventually leads to a very 
efficient reheating of the universe. This is a win-win sit- 
uation, where black holes can reheat the universe even if 
the standard reheating mechanism is inoperative! 

VII. "NATURAL" HYBRID INFLATION 

As we have seen, the model ([[]), as well as the version 
proposed in jL8| , lead to a copious formation of huge black 
holes if one requires that (unlike in the original version of 
hybrid inflation) there is an additional inflationary stage 
after the phase transition. This problem occurs because 
typical classical trajectories in this model go very close to 
-0 = 0. One can avoid this problem by a modification of 
the shape of the effective potential [^3| . Also, as we have 
shown above, black hole production can be even useful if 
the second inflationary stage is very short and the black 
holes are very small. But there exists another problem, 
which we will consider now together with the first one. 

The main reason why many authors are trying to im- 
plement hybrid inflation in supersymmetric theories is to 
protect the flatness of the effective potential in the in- 
direction. One may try to relate the small mass of the 
field <fi to the gravitino mass 7713/2 ~ 1 TeV, which ap- 
pears because of supersymmctry breaking. If one argues 
that the parameter M 2 in Eq. (ph is of the order of the 
intermediate scale of supersymmetry breaking 1113/2 Mp, 
then there appears to be no unexplained small param- 
eters in the model. Still the appearance of the term 
M 4 ~ (m 3 / 2 Mp) 2 cos 2 (i/'/v / 2/) in (|j) remains somewhat 
unclear to us. If one expands M 4 cos 2 (V'/v / 2/) in pow- 
ers of V for M ~ 10 11 GeV and / - 10 18 GeV as in 
Rcf. |18| ], one would find an extremely small coupling 
constant M 4 // 4 ~ 10 -30 in front of the term ip A . It was 
pointed out in | fl8| that such couplings may appear in 
a natural way if one introduces certain superpotentials 
which lead to nonrenormalizable interactions. However, 
to study nonrenormalizable terms in an internally consis- 
tent way it would be necessary to consider models based 
on supergravity, which was outside the scope of our inves- 
tigation, as well as of the investigation performed in |jl8| . 

Fortunately, both smallncss of the parameter M 4 and 
the shape of the potential can be explained if one in- 
terprets ip as a. pseudo Goldstone field similar to the 
axion field. One may consider a model of a complex 
scalar field ^(x) = (f(x)/y/2) expi9(x), which after 
spontaneous symmetry breaking can be represented as 
(f/y/2) exp(iip(x)/f). If the original effective poten- 
tial was a function of \&*\E', the field ip will be mass- 
less. However, nonperturbative (instanton or wormholc) 
effects may give this field a small mass (see Ref. Eq] 
for a recent discussion of this issue). This effect can 
be described by adding operators breaking initial U(l) 
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symmetry of the effective potential. Consider the fam- 
ily of operators g n (*& ± fl/*)™^ -71 . Since these operators 
appear because of nonperturbative effects, the coupling 
constants g n may be exponentially small. One may take, 
for example, the simplest operator gi(^±^>*)f 3 , and add 
to it a constant term y/2g\f A normalizing the vacuum en- 
ergy to zero. This gives the effective potential of the field 
ip, which is completely analogous to the standard axion 
potential,^ 



A! 



V{il>) = 2v / 25i/ 4 cos 2 



(71) 



Note that in this potential ip/f is an angular variable 
from to 2ir. This potential coincides with the effec- 
tive potential of the field ip in Eq. (||) up to an obvious 
change 2 v / 2#i/ 4 — > M 4 , / — > /V2. (Our definition of 
/ corresponds to a canonical normalization of the field 
cp kinetic terms.) In this context both the shape of the 
potential for the field ip and the smallness of the term 
M A cos 2 {ipj\[2f) are explained in a natural way. Poten- 
tials of this type have been used in "natural inflation" 
3|,p6[. The problem with "natural inflation" 



models |25|,| 

is that for the "natural" value of symmetry breaking 
/ < Mp inflation is too short and the spectrum index 
n is significantly less than 1 . There is no such problem in 
our model; the main purpose of the introduction of the 
field ip is to support inflation before the phase transition 
rather than after it. 

Thus it makes a lot of sense to explore cosine poten- 
tials such as Eq. ( |7l| ) in the context of hybrid inflation. 
But with the interpretation of the field ip as a pseudo 
Goldstone particle, one cannot couple it to the field <p 
in the way proposed in Eq. (||). Now ip is the angu- 
lar part of the field 'J, and one cannot write any su- 
perpotentials for the fields 'J and <p which would result 
in the simple interaction terms ~ ip 2 <p 2 . However, since 
we already reinterpreted the cosine term in Eq. (0) as 
appearing from the anomalous term gi(^f + V?*)/ , we 
can go further and introduce an anomalous interaction 



term 52(^6 



^*e ie ) 2 f 2 4> . Note that we have intro- 



duced here for generality the phase shift 9 between the 
two anomalous terms. Indeed, the origin of the first and 
of the second anomalous terms may be different, and a 
priori one should not expect U(l) symmetry to be bro- 
ken by these two terms in the same way. In what follows 
we will assume that 9 is small. The resulting effective 
potential including the mass term of the field <p can be 
represented in the following form: 



"Note that in [[12j a model of hybrid inflation was studied 
based on a supersymmetric Peccei-Quinn model. In that case 
it was the self-interaction energy of the radial degree of free- 
dom f(x) that gave the constant potential energy in the false 
vacuum. 
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Let us now analyze the shape of this potential and its 
relation to the more usual hybrid inflation potential (Q). 
Consider first the case 9 = 0. At large <p the dominant 
term involving ip in Eq. ( |72] ) is the second one, which 
implies that at large <p the field ip will settle in one of the 
minima at ip/ f = rm, where n is some integer. For odd 
values of n the energy density due to the self-interaction 
cosine-squared term also vanishes and we are left with 
conventional chaotic inflation with V = m 2 cp 2 /2. How- 
ever for even values of n, the self-interaction term is non- 
zero and ip is trapped in a false vacuum, like the model 
in Eq. ([[]). Near ip = the potential in Eq. ( |72] ) is given 
by 
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., V 2 + f0V + ^-0 2 . (73) 



One concludes, that about ip = the bare mass squared 
of the field ip is m 2 = —A 2 / 2 , but that the effective 
mass-squared becomes positive for <p > <p c = A1//A2. 
Near <f> c at ip = the energy density is given by 
2A 2 / 2 (/ 2 + m 2 /4A2). The first term dominates, as in 
the usual hybrid inflation, for m <C A2/, and we have 
exactly the same as the model we have analyzed in the 
preceding sections where 
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and the dimensionless parameters introduced in Sec- 
tion H are given by a = (3/167r)m 2 M 2 /A 2 / 4 , j3 = 
(3/16"7r)Mp// 2 and 7 = A 2 ,. The curvature perturba- 
tions produced at (p > (p c are then given by Eq. (^) . For 
small a we have llO|| 
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in order to agree with the COBE normalization . 

There are two regimes one may consider in this theory. 
First of all, one may assume that the coupling constants 
are not extraordinarily small. Then, as was shown in 
Ref. [iTjj , the conditions m H and Sp/p ~ 10 -5 im- 
ply that there was no second stage of inflation in this 
model, i.e. everything is going on as in the first version 
of the hybrid inflation scenario jlQl . One may consider, 
for example, the following parameters: / ~ 10 16 GeV 
(GUT scale), m = 10 10 GeV (intermediate SUSY break- 
ing scale), and Ai ~ A2 ~ 10~ 3 . In this case all condi- 
tions mentioned above will be satisfied. The parameters 
a ~ 10~ 2 and (3 ~ 10 5 so there will be no second stage of 
inflation and no anomalous black hole production. For 
9 = there will be domain wall production. However, 
the domain walls formed in this scenario are unstable 
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because they are always bounded by strings with ^ = 0. 
Moreover, it is sufficient to consider models with a very 
small non-zero value of 9, so that all the evolution will 
go in one direction, and there will be no domain walls or 
other topological defects. One can easily understand this 
if one takes into account that at large <j> the minimum 
of the effective potential with respect to the field ip is at 
ip = Of ^ 0. Thus it is enough to have 9 > 10~ 15 to 
avoid black hole formation in our model. 

As one might expect, the same model for a different 
choice of parameters can reproduce all the results of the 
model proposed in |Q , including the second stage of in- 
flation after the phase transition. This requires to take 
A 2 ~ 10 -30 , which is an extremely small number. How- 
ever, in our case the existence of this small parameter is 
not surprising, because it could appear due to nonper- 
turbative effects. Typically the value of this parameter 
is suppressed by factors such as exp (— 8ir 2 /g 2 ) where g 
is the gauge coupling constant. In some models this sup- 
pression may not be very significant, but in general this 
suppression can easily give numbers much smaller than 
1CP 30 p8| . In particular, in the usual axion theory with 
/ ~ lO^GeV the corresponding constant is of the order 
of 10~ 130 . From this perspective it is more surprising that 
in this model the coupling constant A| is not required to 
be equally small. This disparity can be easily alleviated 
if one does not insist that the masses of both fields as well 
as the Hubble constant at the end of inflation should be 
of the order of m 3 / 2 . 

For 9 — in this model, just like in the model of 
Ref. |18[ ], one obtains inflating topological defects, very 
large density perturbations on the scale corresponding 
to the moment of the phase transition, and catastrophic 
black hole production. However, it is no longer a generic 
property of the model. Remember that the dangerous 
area of the phase space is located very close to <fi = (p c 
and ip = 0. For example, large density perturbations 
are generated only at ip < H ~ 1CT 15 /, for H ~ 10 3 
GeV and / - 10 18 GeV. It is enough to have 9 > 10" 15 
to avoid black hole formation in our model. Thus for 
generic values of 9 inflationary trajectories never come 
close to ip = 0, and the problem disappears. For non- 
integer values of 9/n hybrid inflation can occur along 
ip = 9 + rm for any integer value of n, with the false vac- 
uum energy density equal to 2X 2 f 4 cos 2 (6*/2) for even n 
or 2\f f 4 cos 2 (9/2) for odd values. 

We do not want to pretend that the "natural" hybrid 
inflation model is necessarily very natural. In order to 
study this question one would have to investigate the 
appearance of the anomalous terms in a more detailed 
way, and to analyse the possible effects of adding a more 
standards terms like ty*^(p 2 . Our main purpose was to 
show that the cosine terms with small coefficients can be 
incorporated into the hybrid inflation models, and that 
it is possible to avoid the problem of large density per- 
turbations and black hole production in these models. 
However, as we have argued in the previous section, pri- 
mordial black holes produced after inflation under cer- 



tain conditions may lead to very interesting cosmological 
consequences. On the other hand, there is a simpler way 
to get rid of the black holes; one may simply return to 
the original hybrid inflation scenario without the second 
inflationary stage. 

VIII. CONCLUSIONS 

Inflationary theory was first proposed about 15 years 
ago, and its main principles are by now well understood. 
It is therefore surprising to see that slight modifications 
in simple and natural models may lead to important and 
sometimes absolutely unexpected consequences. For al- 
most fifteen years we knew that inflation exponentially 
dilutes the density of topological defects, but only two 
years ago did we learn that topological defects may in- 
flate themselves |1S( ]. It was also thought that old infla- 
tion did not work, and chaotic inflation always predicted 
that the universe is flat. A year ago, however, it was 
found that the simplest hybrid model where one takes 
old inflation potential for one of the fields and chaotic in- 
flation potential for another one (even if these two fields 
do not interact with each other!) leads to the universe 
consisting of infinitely many separate universes with all 
possible values of Q < 1 |49|| . Now we encountered one 
more surprising fact. For many years it seemed clear that 
inflation erased all pre-existing inhomogencities and did 
not leave much room for the production of primordial 
black holes, which had been the subject of active investi- 
gation in the end of the 70's. Now we see that in a very 
simple inflationary model one can easily obtain a large 
amount of black holes. They are formed only in a spe- 
cific mass range, determined by the duration of inflation 
after the phase transition. Typically they are huge, but 
depending on the parameters of the model they can be 
very small as well. 

Note that black holes are not necessarily a curse but 
could also be a blessing. In the simplest model stud- 
ied here the probability of black hole formation is not 
suppressed at all, and their number appears to be too 
large. We propose some modifications of this model 
where black hole formation is strongly suppressed. It is 
possible for certain values of the parameters of the model 
to have the right amount of relatively large black holes, 
M ~ 10 15 — 10 30 g, that have not yet evaporated and 
may be responsible for the dark matter in the halos of 
galaxies. In a particular model considered in Sect. 6 we 
have shown that the halo of our galaxy may consist of 
black holes of mass ~ 10 21 g. However, numerical val- 
ues of the masses and abundances of the black holes are 
strongly model dependent. Depending on the parame- 
ters of the models there may be just enough black holes 
to have f2o = 1 in the universe. If the black holes are su- 
permassive, one could speculate about their relation to 
the black holes in the centers of galaxies. 

Suppressing the number of large black holes down to 
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a desirable level requires a certain degree of fine tun- 
ing. But it is relatively easy to make the black holes 
very small and harmless by making the second stage of 
inflation short and by ending inflation at large H . For 
H ~ 10 3 GeV the smallest black hole masses ~ Mp/H 
are still very large, about 10 11 g, and they evaporate very 
late, at t ~ 10 3 s . But if, e.g., one takes the models with 
H ~ 10 14 GeV, one can obtain black holes with a mass 
~ 6 x 10 4 g, evaporating at t ~ 10~ 16 s. Such black 
holes would dominate the universe after their formation 
until evaporation. Evaporation of black holes may lead 
to baryon asymmetry production. During the last fifteen 
years this mechanism of baryogenesis was largely ignored 
since it seemed impossible to produce many small black 
holes after inflation (see however Ref . J2^] ) . We may now 
return to the investigation of this interesting possibility. 
Independently of the issue of baryogenesis, one should 
emphasize that the possibility of the black hole dom- 
inance at the intermediate post-inflationary stage may 
change completely the mechanism of reheating after in- 
flation, which would proceed via black hole evaporation. 

In the course of our work we have further developed 
a method of investigation of density perturbations which 
can be applied even for complicated systems of several 
coupled scalar fields [§9|,E58j. This method is rather sim- 
ple and powerful. It gives analytical results in those cases 
in which the motion in field space is integrable, like in hy- 
brid inflation and in theories with coupled inflaton and 
dilaton fields, but it can be used in a more general con- 
text as well. The method consists of three main parts. 
First of all, for any point in the (0, ip) space one finds 
(either analytically or numerically) an inflationary tra- 
jectory going from this point, and calculates the number 
of e-folds N((f>, ip) for this trajectory. This problem is 
easy to solve numerically even for very complicated po- 
tentials. Then one perturbs the position of the initial 
point (</>, ip) by adding to it inflationary jumps, which 
typically are of the order H/2ir, but may be greater or 
smaller, see Section 3. This gives us the perturbation 
of the number of e-folds SN, which is directly related to 
the density perturbations: Sp/p = (4/9)SN((p, ip) at re- 
entry during the radiation dominated era. Note that the 
resulting density perturbations for a given N (i.e. for 
a given wavelength) will depend on the place (</>, ip) our 
trajectory came from. Thus the remaining step is to eval- 
uate the probability that for a given number of e-folds N 
from the end of inflation the field was at any particular 
point (cp, ip). This problem can be solved by using the 
stochastic approach to inflation M] . This approach tells 
us what is the probability to find density perturbations 
of a given amplitude with a given wavelength]^] Usually 



at the last stages of inflation we have a single inflationary 
trajectory, independent of initial conditions. In our case, 
however, inflationary trajectory is unstable at <p < cp c , 
ip = 0; all trajectories bifurcate due to quantum fluctu- 
ations. Therefore the calculation of the probability dis- 
tribution was necessary to show that a typical amplitude 
of density perturbations produced near the point of the 
phase transition is very large. 

As an important by-product of our investigation we 
have found a new type of inflating topological defect. 
They appear in the models where the curvature of the 
effective potential is somewhat greater than H 2 , but nev- 
ertheless the time necessary for symmetry breaking to 
occur is much greater than H~ l . These defects do not 
inflate eternally and do not form a fractal structure found 
in . Still inflation in the cores of these defects contin- 
ues for a while even after it ends outside of them. As a 
result, they lead to large density perturbations of a spe- 
cific type. Until now inflating topological defects could 
be considered as an interesting but somewhat esoteric 
feature of certain inflationary models. Typically the dis- 
tance from us to these defects was many orders of mag- 
nitude greater than the size of the observable part of 
the universe. They were important for understanding of 
the global structure of the universe, but not of our lo- 
cal neighborhood. However, in hybrid inflation models 
with two stages of inflation these defects are abundantly 
produced at the moment of the phase transition, and 
populate the part of the universe which is accessible to 
our observations. We believe that the new type of inflat- 
ing topological defects deserves separate investigation. It 
would be very interesting to understand whether they 
lead to black hole formation and to explore other possi- 
ble observational consequences of these exotic objects. It 
is amazing that very simple models of two scalar fields 
can exhibit such a rich and interesting behavior! 
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FIGURE CAPTIONS 

Figure 1: The amplitude of quantum fluctuations of 
the symmetry breaking field Sip in units of H/2TT, as a 
function of scale k, evaluated when this scale left the 
horizon. Here k c corresponds to the scale that left the 
horizon during the phase transition, when <p — <p c . The 
figure corresponds to generic values of the parameters, 
(a = .3,/3 = 8). The asymptotic value as k — > oo is 
Coo = 8.91. 



Figure 2: The thick dashed-dotted line corresponds to 
the SN = 1 region within which it is not possible to define 
equal-time hypersurfaces, and density perturbations are 
of order one. We also show a few equal-number-of-e-folds 
contours in field space (ip/H, (p/(p c ), for generic values of 
the parameters, (a — .3,/3 = 8). The dashed sections 
enter the SN = 1 region where due to quantum fluctua- 
tions we cannot associate a definite number of e-folds to 
the end of inflation. 



Figure 3: Contour lines for the dispersion of the ip dis- 
tribution at the phase transition (dashed line) and one 
e-fold after the phase transition (continuous line), com- 
pared with the size of the SN = 1 region, (ip 2 ) 1 / 2 = 
nH/2ns (with n = 1, . . . , 7 from the bottom up) in pa- 
rameter space (a, (3). In the region below the curves, the 
probability distribution for the ip field cannot avoid the 
dangerous SN — 1 region and large-amplitude perturba- 
tions will be expected at scales associated with the phase 
transition. 
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